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ABSTRACT 

Analog  methods  for  measuring  rotational  shaft  velocity 
and  acceleration  are  well  developed,  yet  present  many  draw- 
backs in  practical  use.   A  device  is  proposed  to  digitally 
compute  estimates  of  the  instantaneous  velocity  and  accelera- 
tion of  a  shaft.   The  method  of  attack  used  in  arriving  at 
a  final  design  proceeds  from  development  of  signal  generatina 
equipment  through  velocity-  and  acceleration-processing  cir- 
cuitry.  Computer  simulation  is  used  to  optimize  the  elec- 
tronic circuitry  to  achieve  the  best  accuracy. 
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I .   INTRODUCTION 

Devices  which  accurately  measure  the  instantaneous  rota- 
tional velocity  and  acceleration  of  a  shaft  have  countless 
engineering  applications.   One  of  these  is  the  use  of  wheel 
acceleration  and  velocity  measurements  as  inputs  to  an  air- 
craft servo  braking  system  whose  computer  may  temporarily 
override  the  pilot's  braking  commands  should  the  plane  skid 
during  landing. 

Analog  methods  for  measuring  shaft  speed  and  accelera- 
tion are  well  developed.   A  tachometer  geared  to  the  shaft 
being  tested  develops  a  dc  voltage  proportional  to  the  shaft 
speed  while  a  suitable  resistor-capacitor  bridge  approx- 
imately differentiates  this  voltage  to  yield  an  instantane- 
ous acceleration  voltage. 

The  use  of  analog  techniques  for  shaft  measurements  has 
many  drawbacks;  some  of  the  more  objectionable  ones  are: 

1^   The  number  of  significant  figures  in  measurement 
readings  is  limited  to  the  quality  of  analog  equipment  used 
to  generate  and  display  information. 

2)  Additional  equipment  (such  as  a  tachometer)  added 
to  a  shaft  increases  its  inertia. 

3)  Analog  velocity-measurement  devices  are  often  non- 
linear over  the  useful  range  of  velocity  to  be  measured. 

4)  Mechanical  wear  limits  the  useful  life  and  reli- 
ability of  all  mechanical/electrical  generating  units. 


Since  analog  methods  offer  many  problems,  a  need  was  seen 
for  the  design  of  a  system  to  digitally  measure  the  instanta- 
neous velocity  and  acceleration  of  a  shaft.   This  thesis  is 
the  result  of  the  study,  design,  and  simulation  of  digital 
electronic  equipment  to  achieve  that  end. 

II.   DESIGN  PROCEDURE 

The  design  of  a  digital  system  to  measure  shaft  accel- 
eration and  velocity  was  in  effect  the  design  of  a  small, 
special-purpose  computer  to  perform  these  functions.   The 
design  procedure  for  this  study  was  broken  into  seven 
categories : 

1)  Analysis  of  design  considerations  related  to  rota- 
ting shaft  measurements;  statement  of  formulae  used 
to  express  rotary  motion  and  expression  of  the  re- 
quirements to  be  met  by  the  proposed  equipment 
design . 

2)  Development  of  an  "inertialess"  signai  generating 
source  whose  output  is  directly  proportional  to 
shaft  speed. 

3)  Development  of  digital  circuitry  to  process  the 
shaft  speed  signal  and  produce  velocity  information. 

4)  Development  of  digital  circuitry  to  process  the 
velocity  information  and  produce  acceleration 
information, 

5)  Computer  simulation  of  the  total  digital  design  to 
yield  simulated  data  input  and  digital  system  out- 
put information. 

6)  Analysis  of  the  computer  simulation  to  optimize 
the  system  accuracy. 

7)  Evaluation  of  the  results. 


III.   DESIGN  CONSIDERATIONS 

Necessary  definitions  concerning  the  rotation  of  a 
(rigid)  shaft  can  be  easily  expressed.   In  a  specific  inter- 
val of  time,  a  rotating  shaft  turns  through  a  rotational 
angle,  A0.   The  average  angular  speed,  ou  is  defined 

-    e2-6l 

u>  =  ^i  (1) 

2   1 
where  A0  is  the  rotational  displacement  the  shaft  has  under- 
gone during  time  At.   This  can  also  be  expressed  as 

-     Vei 

w  =  j.  _.,.    •  (2) 

z2    ^1 

If  At  approaches   0   as  a  limit,  instantaneous  angular  speed 

can  be  expressed  as 

_  ,     A9    d9  ,_. 

CO  =  llltl  -~r    =  tt   .  (3) 

At-0  At    dt 
Angular  acceleration,  a,  is  the  time  rate  of  change  of 
angular  velocity.   Average  acceleration,  a,    can  be  expressed 

Au)     2   1  ,  „  v 

where  go-  is  final  and  go,  is  initial  angular  speed  and  t„ 
is  final  time  and  t,  is  initial  time. 

When  dealing  with  problems  concerning  uniformly  accel- 
erated rotary  motion  the  following  equations  apply: 

0  =  cot  +  0.  (5) 

co~  =  co,  +  at  (6) 

02  =  u^t  +  j   at2  +  B1  (7) 


where  t  is  elapsed  time  from  initial  conditions,  and  the 
other  variables  along  with  their  subscripts  are  defined  above. 

The  design  of  a  digital  device  requires  that  the  output 
be  quantized,  and  that  it  appear  at  discrete  output  times. 
The  design  of  this  special-purpose  computer  is  predicated 
upon  the  successful  implementation  of  electronic  devices  to 
perform  calculations  in  equations  (2)  and  (4) 

-    02-0l 

a)  =  ^i  (2) 

r2  rl 

00^-00, 

a  =  r^rr^   •  (4) 

c2  ^1 

Although  it  would  have  been  desirable  to  design  a  system 
able  to  measure  any  range  of  rotational  speeds,  it  was  felt 
that  design  considerations  should  be  predicated  on  solving  a 
specific  problem.   The  specific  problem  undertaken  was  the 
design  of  a  special-purpose  digital  computer  to  measure  the 
rotational  speed  and  acceleration  of  a  magnetic  tape  trans- 
port capstan.   The  tape  transport  system  concurrently  under 
design  by  a  computer  company  was  to  provide  an  unusual 
acceleration:   0   through  200  inch/second  tape  speed  over 
elapsed  time  of  3.75  msec  using  a  capstan  of  0.75  inch 
radius.   To  insure  that  this  thesis  design  might  adequately 
handle  this  problem,  the  tape  drive  requirements  were  made 
more  demanding. 
Tape  Drive  Requirements: 

Drive  capstan  radius  0.75  inch 

Operational  tape  speed  200  inches/sec 

Acceleration  time  (0-200  inch/sec)       2.5  msec 
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These  data  were  used  to  compute  average  shaft  speed  and 

acceleration.   The  average  linear  velocity,  v,  of  the  capstan 

surface  can  be  expressed 

v  =  rco  (8) 

where  r  is  the  capstan  radius  and  co  is  average  shaft  speed. 

Converting  this  equation  and  substituting  values  yields 

v    200  inch/sec  x  1  revolution    „ »  _  rev   /riN 

(jo  =  —  =  — = — =-= - — r ^ 3-1 =  42.5  .  (9) 

r      0.75  inch  x  2 tt  radians  sec 

The  average  acceleration  can  be  computed  by  equation  (3). 

2   1    42.5  rev/sec  -  0  rev/sec    .  _    , „ 4  rev 
2   1     2.5  x  10    sec  -  0  sec  sec 

(10) 

Allowing  for  some  overdesign  it  was  required  that  the 

system  be  able  to  measure  (at  least)  the  following  parameters 

Expected  Shaft  Speed  and  Acceleration: 

4         2 
max  acceleration  -  2  x  10   rev/sec 

max  speed  -  100  rev/sec  . 

IV.   DEVELOPMENT  OF  A  SIGNAL  GENERATING  SOURCE 

The  many  inertialess  signal  generating  methods  studied 
were  capable  of  generating  electrical  signals  directly  pro- 
portional to  shaft  position  or  velocity-   A  photoelectric 
design  most  suitable  to  the  problem  solution  is  hereafter 
described- 
It  was  anticipated  that  the  capstan  drive  shaft  could 
have  black  lines  (arc  of  0.2°  width)  easily  scribed  every 
3.6  degrees  around  its  circumference.   One  small  photocell 
measuring  the  intensity  of  reflected  light  from  these  lines 
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would  produce  a  resistance  inversely  proportional  to  the 
square  root  of  that  intensity.   To  increase  the  resolution  of 
this  system  (from  3.6°  to  0.36°)  the  following  scheme  was 
proposed. 

As  shown  in  Figs.  1  and  2,  there  are  lines  scribed  on  the 
capstan  above  which  is  positioned  a  photocell  mountinq  brack- 
et.  Nine  photocell  "windows"  (0.2  degrees  of  arc  wide)  are 
cut  in  this  bracket  at  an  equal  radius  from  the  capstan  cen- 
ter and  equally  spaced  precisely  four  degrees  apart. 

Looking  at  Fig.  2  one  can  see  that  only  one  photocell 
window  has  a  line  completely  under  it  at  any  time.   If  the 
capstan  is  turned  clockwise  0.36  degrees  it  is  easily  seen 
that  photocell  number  two  would  have  directly  under  it  the 
line  which  previously  had  been  under  photocell  number  one. 
The  resolution  of  the  system  is  therefore  0.36°.   A  shaft 
rotation  of  3.6°  would  cause  ten  sequential  signals  to  be 
developed  due  to  the  vernier  construction  of  this  device. 


The  photocell  envisaged  for  this  system  is  listed  in 
Industrial  Catalog,  number  680,  p.  83,  Allied  Electronics, 
1968.   An  "NPN  Planar  Silicon  Light  Sensor  #  LS-400"  has  the 
following  response  times:  enable  time  1.5  psec;  disable  time 
15  usee.   Assuming  a  total  cycle  time  of  20  ysec,  the  maximum 
shaft  rotational  speed  with  errorless  signal  generation  can 
be  easily  calculated  for  this  system. 

-. r-TTTr — =-: = r -,  .-,  ■« >  20  ysec 

x  rev/sec  •  100  lines/rev  •  1  output  signal/line  — 

(11) 

Solving  for  x  yields 

X  ±  100  .  20  *  10-*  =  50°  rev/sec         <12> 

This  photocell  would  easily  fulfill  the  requirements  of  the 
tape  drive  system. 
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FIG.  1   ENLARGED  VIEW  OF  SIGNAL  GENERATOR 
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(TOP  VIEW) 
FIG.  2   ENLARGED  VIEW  OF  SIGNAL  GENERATOR 
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The  photocell  of  a  "Gossen,  Lunasix"  photo  lightmeter  was 
independently  studied  and  found  to  vary  between  70  ohms  and 
15  kilohms  depending  on  the  amount  of  light  impinging  on  the 
cell.   Photocells  for  the  system  under  design  could  be  ex- 
pected to  act  similarly.   It  was  anticipated  that  all  nine 
photocells  could  be  placed  in  parallel  and  input  to  the  cir- 
cuit shown  on  the  following  page. 

A  description  of  the  electronic  operation  follows.   An 
incremental  (0.36  degrees)  shaft  displacement  signal  is  devel- 
oped across  resistor,  RIO  and  the  photocells  in  parallel. 
Coupled  by  CI  to  the  base  of  Ql,  the  positive  bias  of  Ql  is 
overcome  and  the  transistor  is  temporarily  cut  off.   The 
amplified  signal  is  coupled  by  C2  to  a  multivibrator  composed 
of  the  remaining  elements.   R15  serves  as  a  pulsewidth  ad- 
justment, whose  time  constant  is  determined  by  the  product  of 
R15  x  C3.   Diode  D2  serves  as  a  voltage  stabilizing  diode. 

In  summation,  circuitry  has  been  developed  to  generate  a 
variable  pulse  width,  positive-going  square  wave  that  is 
generated  whenever  the  shaft  is  rotated  every  0.36  degrees. 

V.   DEVELOPMENT  OF  CIRCUITRY  TO  PROCESS  THE  SHAFT 
SIGNALS  AND  OUTPUT  VELOCITY  INFORMATION. 


The  development  of  digital  circuitry  to  process  the  shaft 
positional  signals  involved  three  important  concepts. 

1)  The  shaft  signal  had  to  be  stored  or  "buffered"  until 
the  processing  circuitry  had  processed  it. 
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2)  An  electronic  method  of  digitally  generating  the 
velocity  signal  had  to  be  developed. 

3)  After  the  new  velocity  signal  had  been  generated, 
circuitry  had  to  be  designed  to  remove  the  previous  signal 
stored  in  the  buffer  to  enable  it  to  receive  another  outout 
from  the  signal  generator. 

The  final  circuitry  anticipated  to  fulfill  the  require- 
ments above  is  shown  on  the  next  page. 

The  digital  components  proposed  for  this  accelerometer 
are  basic  logic  elements.   The  design,  theory  of  operation, 
and  Boolean  algebra  representation  of  the  basic  elements 
used,  such  as  AND  gates,  OR  gates,  JK  flip-flops,  and  shift 
registers  are  described  in  many  textbooks  [Refs.  1,  3,  and 
5],   Other  equipment  such  as  clock  generators,  divide-by-500 
counters,  decoders,  and  output  display  equipment  are  fairly 
standard.   Since  the  proposed  accelerometer  was  not  con- 
structed, a  detailed  circuit  design  has  not  been  included. 

An  explanation  regarding  individual  component  functions 
in  relation  to  system  operation  can  now  be  made.   During  the 
operation  of  the  accelerometer  there  is  received  a  non- 
periodic,  positive-going  square-wave  signal  for  each  incre- 
mental shaft  rotation  of  0.36°.   If  the  normal  output  state 
for  the  JK  flip-flop  is  a  "0",  the  occurrence  of  the  square- 
wave  signal  at  the  set,  J  terminal  causes  a  "1"  to  occur  at 
the  output  terminal,  Q.   (This  "1"  output  remains  until  the 
K  terminal  has  a  "1"  input  into  it) . 
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After  the  flip-flop  output  is  "1"  the  occurrence  of  a 
positive  clock  signal  causes  the  output  of  AND  gate  number 
one  to  be  "1".   The  use  of  AND  gates  numbered  two  and  three 
permits  a  sufficient  time  delay  to  occur  such  that  AND  gate 
number  one  has  reacted  long  enough  to  guarantee  output  of  a 
"1"  signal.   This  delayed  clock  enables  the  1-bit  shift  reg- 
ister (which  already  has  a  "1"  input  to  the  J  input  terminal) 
and  the  output  shifts  from  an  "0"  to  a  "1".   The  delayed 
500-kHz  clock  is  also  fed  to  AND  gate  number  four.   When  this 
delayed  clock  and  the  output  of  the  1-bit  shift  register  are 
both  "1"  AND  gate  number  four  outputs  a  "1",  which  resets  the 
JK  flip-flop  to  the  "0"  output  state.   --  An  important  design 
requirement  should  be  noted.   The  output  from  the  shaft  sig- 
nal generator  must  be  less  than  1  microsecond  in  duration 
(half  the  clock  period)  to  eliminate  the  chance  of  propa- 
gating two  output  signals  to  the  shift  register  for  one  in- 
put shaft  signal.   --  Thus  circuitry  has  been  designed  to 
buffer  the  shaft  signal  until  it  is  processed  and  sent 
through  the  1-bit  shift  register,  at  which  time  the  signal 
is  removed  from  the  flip-flop,  preparing  it  to  receive 
another  pulse. 

From  the  circuitry  described  it  can  be  seen  that  the 
1-bit  shift  register  has  0's  or  l's  propagated  through  it 
at  every  clock  time,  a  period  of  two  yseconds.   From  Fig.  5 
it  can  be  seen  that  the  1-bit  shift  register  is  followed  by 
a   100-bit  shift  register.   The  time  for  a  single  "1"  to 
propagate  through  the  long  shift  register  can  be  calculated. 


Propagation  Time  =  number  of  shift  register 

bits/clock  frequency         (13) 

Thus  if  a  200-bit  shift  register  were  used  with  a  clock  fre- 
quency of  500  kHz 

Propagation  Time  =  200/5  x  10   =  400  psec  .    (14) 
The  total  number  of  l's  in  the  shift  register  is  indic- 
ative of  the  average  velocity  sensed  over  the  propagation  time 
interval.   An  example  best  demonstrates  this. 

Assume  the  number  of  shift  register  bits  is  500,  and  that 
the  clock  frequency  is  500  kHz.   Assume  also  that  the  signal 
sensor  has  been  sensing  signals  from  a  shaft  which  has  been 
continuously  rotating  at  10  revolutions/second.   The  follow- 
ing calculations  can  be  made. 

10  revolutions/second  x  360  degrees/revolution  x  1 

signal/0.36  degrees  =  10,000  signals/second.    (15) 

Propagation       500  shift  register  bits _  -,  nnn 

5  x  10   cycles/sec  x  1  shift/cycle 

(16) 

4  -3 

10   signals/second  x  10    sec  =10    l's  stored  in  the 

shift  register  . 

(17) 
From  this  example,  the  number  of  l's  stored  in  the  loner  shift 
register  directly  corresponds  to  the  shaft  velocity  being 
measured.   The  method  of  adding  the  l's  in  the  shift  register 
is  effectively  handled  by  the  use  of  an  "up-down"  counter, 
a  device  which  has  the  capability  of  increasing  or  decreas- 
ing its  count  when  the  appropriate  "up"  or  "down"  command  is 
present.   The  occurrence  of  the  delayed  500  kHz  clock  causes 
the  serial  shift  register  to  simultaneously  shift  right  by 
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one  bit  position  all  of  its  stored  bits.  AND  gates  five  and 
six  are  included  to  prevent  "up"  and  "down"  signals  arriving 
simultaneously,  a  condition  not  permitted  in  up-down  counters 

Circuitry  developed  to  this  point  has  stored  in  the  up- 
down  counter  the  total  number  of  electronic  signals  corre- 
sponding to  an  average  shaft  velocity.   These  signals  are 
sensed,  decoded  (to  decimal  numbers),  and  displayed  at  every 
clock  time.   Thus  circuitry  has  been  designed  to  process  the 
shaft  speed  signal  and  display  velocity  information. 

VI.   DEVELOPMENT  OF  CIRCUITRY  TO  PROCESS 
ACCELERATION  INFORMATION 

From  equation  (4)  it  can  be  seen  that  average  accelera- 
tion over  an  interval  is  the  change  in  angular  shaft  veloc- 
ity divided  by  the  time  interval.   The  development  of 
digital  circuitry  to  perform  this  calculation  involved 
three  important  concepts. 

1)  An  electronic  method  of  generating  the  acceleration 
output  period  had  to  be  developed. 

2)  The  shaft  velocity  had  to  be  sampled  and  stored 
once  each  period. 

3)  An  electronic  method  for  calculating  and  displaying 
the  acceleration  value  once  each  period  had  to  be  devised. 

Fig.  4  includes  the  circuitry  to  perform  these  functions. 
Individual  component  operation  can  now  be  discussed. 

Assume  a  shaft  velocity  has  been  calculated  and  is  there- 
fore stored  in  the  up-down  counter.   Assume  also  that  at  time 
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"tl"  a  velocity  has  been  stored  in  both  8-bit  shift  registers 
The  divide-by-500  counter  effectively  generates  a  1000-Hz 
clock  frequency  that  is  synchronized  with  the  master  (500kHz) 
clock.   When  the  next  1000  Hz  clock  pulse  is  generated  at 
time  "t2",  this  velocity  is  shifted  (after  AND  gate  number 
nine's  time  delay)  into  the  positive  terminal  of  the  sub- 
tractor.   After  AND  gate  number  seven's  time  delay  the  new 
velocity  in  the  up-down  counter  is  sampled  and  stored  in  the 
8-bit  shift  register.   After  AND  gate  number  eight's  time 
delay  this  velocity  is  shifted  simultaneously  to  the  output 
terminals  of  the  8-bit  shift  register  one  as  well  as  to  the 
minus  terminals  of  the  subtractor.   After  AND  gate  number 
10 's  time  delay  the  subtractor  is  enabled,  permitting  the 
subtraction  of  velocities  sampled  at  times  "t2"  and  "tl". 
The  input  to  the  acceleration  decoder  is  therefore  the  dif- 
ference in  velocities  sampled  at  millisecond  intervals.   The 
value  1000  as  a  scale  factor  must  be  used  to  yield  true 
acceleration. 

Circuitry  has  been  devised  to  calculate  average  shaft 
acceleration  at  1-millisecond  intervals.   Although  informa- 
tion displayed  at  this  rate  is  too  fast  for  human  observa- 
tion, it  is  anticipated  that  this  acceleration  information 
would  be  electronically  incorporated  with  a  servo  control 
system. 
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VII.   COMPUTER  SIMULATION  OF  THE  SYSTEM 


It  was  anticipated  that  the  system's  calculations  of 
velocity  and  acceleration  would  have  some  error  due  to  the 
quantizing  necessary  in  a  digital  system.   It  was  expected 
that  these  errors  would  have  some  complicated  relationship 
to  the  number  of  shift  register  bits  used  to  store  the  l's 
and  O's  corresponding  to  shaft  velocity.   A  computer  simula- 
tion was  written  to  experimentally  determine  this 
relationship. 

A  printout  of  the  simulation  used  is  printed  under  the 
Computer  Simulation  Program  Listing  section.   The  most  impor- 
tant simulation  variable  names  are  listed  below  with  a 
description  of  their  function. 


Variable 


TIME 


DELT 


FINTIM 


CKOUT 


THETA 


OMEGA 


Purpose  Units 

Serves  as  the  independent  time 

base  for  the  simulation.  seconds 

Serves  as  the  incremental  time 

base  for  the  simulation.  seconds 

Serves  as  the  final  time  to  be 

simulated.   The  total  real  time 

simulation  was  5.15  milliseconds.     seconds 

Serves  as  the  master  500kHz  clock; 
it  is  a  square  wave  of  amplitude 
1.0  and  has  a  period  of  two 
microseconds . 

Serves  as  the  true  angular  position 
of  the  simulated  shaft  which  is 
uniformly  accelerated  from  rest  in 
accordance  with  equation  (7).         degrees 

Serves  as  the  instantaneous  true 

shaft  velocity.  degrees/second 
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Variable 


Purpose 


Units 


OMi-GAO     Serves  as  the  initial  shaft 
velocity  at  time  zero. 

ALFA       Serves  as  the  true  simulated 
acceleration. 


degrees/second 


degrees/second' 


ALRPS2     Serves  as  the  true  simulated 
acceleration . 


r evolution s/ second' 


RPST 


Serves  as  the  true  velocity 
of  the  simulated  shaft. 


revolutions/second 


SPIKE      Serves  as  the  electrical  out- 
put voltage  developed  across 
the  nine  photocells  in  the 
Signal  Generator.   -There  is 
a  one- volt  "spike"  developed 
when  the  shaft  has  turned 
each  0.36  degrees, 

SHIN      Serves  as  the  electrical  volt- 
age output  as  a  0 . 9-millisecond 
positive-going  square  wave  from 
the  Shaft  Position  Generator 
Circuitry  shown  in  Fog.  3. 

FFOUT      Serves  as  the  binary  output  from 
the  JK  flip-flop.   fFFOUT  is 
either  0  or  1 . ) 

IDD2       Serves  as  the  delay  time  required 
for  the  JK  flip-flop  to  change 
its  output  from  0  to  1 .   (IDD2 
expresses  in  integer  form,  the 
number  of  DELT  simulation  times 
contained  in  an  average  JK  flip- 
flop  enable  time.) 

IDD3      Serves  as  the  delay  time  required 
for  a  JK  flip-flop  to  change  its 
output  from  1  to  0 .   (IDD3  ex- 
pressed in  integer  form,  the  num- 
ber of  DELT  simulation  times 
contained  in  an  average  JK  flip- 
flop  disable  time.) 

BFOUT      Serves  as  the  output  from  AND 
gate  number  1 . 

IDD4       Expresses  the  integer  number  of 

DELT  simulation  times  contained  in 
an  average  AND  gate  enable  time. 


volts 


volts 


2.3 


Variable  Purpose 

IDD5       Expresses  the  integer  number 
of  DELT  simulation  times  con- 
tained in  an  average  AND  gate 
disable  time. 


Units 


DCK       Serves  as  the  delayed  500kHz 
clock  after  going  through  AND 
gates  two  and  three. 

SR10UT     Serves  as  the  output  from  the 
1-bit  shift  register. 


SHRSET 


NUM 


Serves  as  the  output  from  AND 
gate  number  four;  its  effect 
is  to  reset  the  JK  flip-flop 
to  output  0  and  enable  it  to 
receive  another  SHIN  signal. 

Serves  as  the  integer  number 
of  shift  register  bits  being 
simulated. 


MINUS      Serves  as  the  output  from  the 
100-bit  shift  register.   When 
MINUS  is  1,  the  total  stored 
in  the  up-down  counter  is  de- 
creased by  a  value  of  1. 


RPSE 


Serves  as  the  shaft  velocity 
as  stored  in  the  up-down 
counter. 


revolutions/second 


STORE1     Serves  as  the  velocity  stored 
in  the  8-bit  shift  register 
number  1 . 


re volutions /second 


STORE2     Serves  as  the  velocity  stored 

in  the  8-bit  shift  register    revolutions/second 
number  2 . 


ACCEL      Serves  as  the  output  from  the 
acceleration  decoder. 


re volutions/ second 


MAINN      A  subroutine  which  serves  as 
the  main  program.   The  mam 
program  (which  calls  MAINN) 
serves  as  a  dummy  program  whose 
purpose  is  to  dimension  a  matrix 
which  corresponds  to  the  100-bit 
shift  register. 
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Variable  Purpose  Units 

STPLOT     A  subroutine  which  graphically 
displays  system  responses. 
(The  sample  time  used  for  this 
simulation  is  1  sample/20 
microseconds  to  permit  photo- 
graphic reduction  of  the  re- 
sults for  inclusion  in  the 
thesis . ) 


VI.   EVALUATION  OF  THE  RESULTS 

Fifteen  computer  calculations  were  made  to  simulate  the 
system  output  while  using  the  following  number  of  shift 
register  bits:   32,  40,  64,  84,  100,  108,  120,  130,  150, 
200,  250,  350,  400,  500,  and  640.   Average  velocity  disper- 
sion data  were  derived  by  drawing  two  straight  lines  through 
the  most  dispersed  "RPSE"  points  on  the  computer  graphs,  Fig. 
7.   The  vertical  distance  between  these  lines  was  averaged 
(for  1,  2,  3,  4,  and  5  millisecond  times)  and  divided  by 
two  to  yield  the  average  maximum  distance  from  the  mean. 
Two  straight  lines  were  drawn  through  the  top  and  bottom 
points  of  the  character  used  to  plot  the  curve  "RPST".   The 
vertical  distance  between  these  two  lines  was  subtracted 
from  the  average  maximum  distance  derived  above.   This  re- 
maining distance,  when  multiplied  by  the  proper  graph  scale, 
10.0  units/inch,  yielded  the  average  velocity  dispersion. 
Two  lines  were  plotted  midway  between  each  pair  of  dispersion 
lines  mentioned  above  yielding  accurate  graphical  depictions 
of  true  and  estimated  velocity.   The  vertical  distance  be- 
tween these  two  lines  was  averaged  and  scaled  to  yield  the 
average  velocity  error.   Average  acceleration  data  were 
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derived  from  averaging  the  written  computer  acceleration  out- 
put values  computed  after  one  millisecond.   The  averaged 
error  data  collected  from  this  analysis  are  listed  below. 


Number 

Average 

Average 

Average 

of  Shift 

Velocity 

Velocity 

Acceleration 

gister  Bits 

Error 
(Rev/Sec) 

Dispersion 
(Rev/Sec) 

2 
(Rev/Sec  ) 

32 

-  3.300 

13.85 

1.875  x  104 

40 

-  3.180 

11.35 

1.750 

64 

-  2.862 

6.26 

1.875 

84 

-  2,940 

5.31 

1.905 

100 

-  2.531 

3.72 

1.900 

108 

-  3.158 

3.58 

1.945 

120 

-  3.280 

3.17 

1.917 

130 

-  3.412 

2.86 

1.923 

150 

-  3.593 

2.30 

1.933 

200 

-  4.891 

1.90 

1.960 

250 

-  4.500 

1.48 

1.880 

350 

-  6.844 

1.08 

1.857 

400 

-  8.688 

0.91 

1.850 

500 

-10.063 

0.53 

1.800 

640 

-11.187 

0.52 

1.750 

These  data  are  depicted  graphically  in  Fig.  5. 

An  evaluation  of  the  computer  results  can  now  be  made. 

Fig.  5  shows  that  the  average  acceleration  peaks  close  to  an 

4  2  4 

average  of  1.9  x  10   revolutions/second   (versus  2  x  10 

which  was  simulated)  for  a  shift  register  bit  length  of  100. 
The  average  velocity  error  reaches  a  minimum  and  has  a  veloc- 
ity dispersion  of  less  than  four  revolutions/second  for  this 
same  number  of  shift-register  bits.   For  these  reasons,  it  is 
recommended  that  the  final  shift-register  bit  length  be  made 
one  hundred. 

Two  amplifying  points  should  be  made  to  stress  this  deter- 
mination.  The  high  dispersion  caused  by  simulation  of  a  40- 
bit  register  caused  a  large  acceleration  error.   An  error  such 
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as  this  would  stress  the  point  of  using  a  low  number  of  shift- 
register  bits.   The  choice  of  using  100  bits,  on  the  other 
hand,  appears  to  yield  a  better  calculation  of  acceleration. 

Although  not  shown  in  Fig.  5,  the  acceleration  calculated  by 

4 
the  system  for  a  100-bit  shift  register  is  exactly  2.0  x  10 

for  the  last  4  milliseconds. 

As  a  final  evaluation  of  the  system  a  frequency  response 

of  the  velocity  computing  circuitry  was  simulated  on  the 

2 
computer.    A  summary  of  the  changes  made  is  listed  below. 

1)  Initial  shaft  velocity  was  42.5  revolutions/second. 

2)  The  shaft  velocity  was  sinusoidally  modulated  about 
this  average  shaft  speed  at  frequencies  of  interest. 

3)  The  instantaneous  accelerations  needed  to  achieve 
these  shaft  velocity  changes  were  calculated  and 
plotted 


2 
It  was  unfeasible  to  make  a  frequency  response  simula- 
tion of  the  acceleration  computing  circuitry  due  to  two 
factors.   Since  the  acceleration  calculations  were  made  for 
every  millisecond  of  simulated  time,  a  significant  number  of 
individual  acceleration  calculations  was  limited  to  the 
amount  of  real  computer  time  available  for  an  extended  study. 
The  amount  of  significant  experimental  evidence  to  be  gained 
from  at  least  70  minutes  of  anticipated  Central  Processing 
Unit  time  on  the  IBM  360/67  was  considered  too  expensive  to 
justify  the  study.   More  important,  however,  is  the  fact  that 
the  acceleration  sampling  rate  imposes  a  limit  on  the  varia- 
tions that  can  be  measured  accurately.   Specifically,  C.  E. 
Shannon's  Sampling  Theorem  states  that  "the  sampling  frequency 
must  be  at  least  twice  as  great  as  the  frequency  of  the  high- 
est Fourier  component  of  the  analog  signal  being  sampled." 
[Ref.  2]  in  other  words,  this  accelerometer  cannot  be  expected 
to  calculate  accurately  accelerations  which  cause  a  shaft  to 
change  its  response  at  a  rate  exceeding  500  Hz. 
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4)   The  input  shaft  velocities  and  the  system  responses 
were  simultaneously  plotted  for  later  frequency 
response  calculations. 

The  actual  computer  simulation  changes  used  are  listed 

below. 

Change  Computer  Format  Card  Number 

1.  DATA  TIME,  OMEGAO , START/0 . ,  15  300.,  0./   *   00027 

2.  FREQMD  -  (the  modulation  frequency        *   0  0  066 

simulated) 

3.  OMEGA1  =  OMEGA  *   022 7 A 

4.  OMEGA  =  OMEGAO  *  COS  (TWOPI *FREQMD*TIME) 

+  OMEGAO  *   00233 

5.  ALFA  =  (OMEGA  -  OMEGAl)  /  DELT  *   02  33A 

6.  2  T79,  ' RPST',  T88,  'ALFA'  /  )  *   00388 

7.  4000  FORMAT (IX,  16,  1PE13.6,  T34 ,  1PE13.6, 

T55,  *   00394 

8.  *  1PE13.6,  T79,  1PE13.6,  2X ,  1PE13.6)     *   0394A 

Ten  computer  simulations  were  made  to  determine  the  re- 
sponse for  the  velocity  computing  circuitry.   The  results  are 
tabulated  below  and  graphed  on  the  next  page. 


Modulation 

Phase 

Frequen 

cy 

Response 

Lag 

(Hz) 

(decibels) 
0. 

(degrees) 

20 

undeterminable 

60 

0. 

less  than  2.0 

150 

0. 

6.42 

500 

0. 

17.9 

850 

0. 

32.2 

1300 

-  0.54 

43.5 

1800 

-  0.54 

63.5 

2000 

-  1.17 

66.9 

2500 

-  1.89 

90.0 

3500 

-  3.85 

119.8 

5000 

-12.24 

undeterminable 
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Fig.  7  demonstrates  well  that  the  system  has  a  usable  out- 
put far  in  excess  of  the  500  Hz.  restriction  imposed  by  the 
acceleration  sampling  rate.   It  should  be  noted  that  the  shaft 

acceleration  necessary  to  produce  a  500  Hz.  modulation  is 

5  2 

1.334  x  10   revolutions/second  .   Since  Fig.  7  shows  that  the 

proposed  system  can  perform  well  for  modulation  frequencies 

of  500  Hz. ,  the  anticipated  use  of  velocity  and  acceleration 

calculations  which  vary  at  much  slower  rates  is  expected  to 

be  reliable. 
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